Screening model for nanowire surface-charge sensors in liquid 
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The conductance change of nanowire field-effect transistors is considered a highly sensitive probe 
for surface charge. However, Debye screening of relevant physiological liquid environments challenge 
device performance due to competing screening from the ionic liquid and nanowire charge carriers. 
We discuss this effect within Thomas-Fermi and Debye-Hiickel theory and derive analytical results 
for cylindrical wires which can be used to estimate the sensitivity of nanowire surface-charge sen- 
sors. We study the interplay between the nanowire radius, the Thomas-Fermi and Debye screening 
lengths, and the length of the functionalization molecules. The analytical results are compared to 
finite-element calculations on a realistic geometry. 



Imagine a sensor so small and compact that it can 
fit almost everywhere and detect all kinds of chemical 
substances in real time. Such a sensor could in princi- 
ple monitor and detect unwanted bacteria and viruses 
instantaneously in eg. your blood or drinking water. 
Sensors based on semiconductor nanowires have already 
been fabricated [l[ and shown to work, for example as 
a pH sensor, where the concentration of hydrogen ions 
H + in a surrounding liquid is detected 0- Moreover, 
applications in label-free detection and biological sensing 
addressing e.g. DNA in low concentration is now being 
explored [U, [j, |5, 6]. In general, the conductance G of a 
nanowire is in the literature considered a promising can- 
didate for a highly sensitive probe of charged particles 
covering or situated near its surface. Detection limits 
have been studied theoretically with respect to binding- 
diffusion dynamics of the molecules at and near the sur- 
face @, However, it is also crucial to address the fun- 
damental problem of effective screening of these charges 
in physiologically relevant liquids which typically yield 
screening of charges on the nanometer length scale. Of 
course it is possible to decrease the salt concentration in 
the analyte solution, but this may not only complicate 
the process - it may also change the biological functions. 

Effects of screening have been discussed qualitatively 
in the literature (see e.g. Refs. [H, touching on the in- 
fluence from the competing screening mechanisms of the 
electron gas in the nanowire and the dissolved ions in 
the surrounding electrolyte. Very recently the effect has 
also been discussed in a computational study of a Silicon 
nanowire (SiNW) FET |Ioj]. In this Letter we develop 
a simple screening model which exemplifies how the de- 
sired screening of surface charge by the nanowire charge 
carriers may be jeopardized by the contending screen- 
ing by the electrolyte. We derive analytical results for 
cylindrical wires which can be used to estimate the sen- 
sitivity of a nanowire surface-charge sensor surrounded 
by electrolyte. The analytical results are compared to 
finite-element calculations on a realistic geometry 

A surface charge density as will perturb the initial 
charge-carrier density no in the nanowire by An, thus 
changing the initial conductance Go by AG. As an ex- 
ample, a positive surface charge will attract additional 



electrons from the nanowire contacts (equilibrium reser- 
voirs). From simple conservation of charge we arrive at 
the following expression for the sensitivity 
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with r being a dimcnsionless function, between zero and 
unity, quantifying the actual sensitivity in the presence of 
Debye screening in the electrolyte and a finite Thomas- 
Fermi screening in the nanowire. We show that T scales 

1/3 

with the electron density no as n in the dilute limit, 

2/3 

thus leading to a prediction of a n dependence, con- 
trasting the intuitively expected n^ 1 dependence. How- 
ever, high sensitivity is of course still associated with low 
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FIG. 1: Schematic cross section of a nanowire covered by 
an ultra-thin oxide (<5) immersed in an electrolyte. A posi- 
tive surface- charge density as, supported by functionalization 
molecules of length £, is screened by negative ions from the 
electrolyte and electrons from the nanowire. For the opposite 
sign of as, the surface charge will be screened by attracting 
positive ions from the electrolyte and by expelling the elec- 
trons from the nearby-surface region of the nanowire. 
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densities. 

Above, P is the perimeter of the of nanowire cross 
section (the fraction supporting the surface charge), A 
is the cross-sectional and e is the electron charge. 

For a cylindrical wire of radius R we have P/A = 2/R, 
thus clearly illustrating the benefit of scaling the wires 
to the nano regime. Obviously, screening in the liquid 
will suppress the sensitivity below the bound given by 
Eq. |T]) and a v,q 1 dependence (the limit T = 1) can 
only be expected when the Thomas-Fermi screening in 
the nanowire is much stronger than the Debye screen- 
ing in the electrolyte surrounding the nanowire so that 
changes in the density of the electron gas fully compen- 
sates the additional surface charge. In our screening 
model we consider the induced electrical potential due to 
a surface-charge density as on the outside of an oxide- 
covered nanowire, see Fig.Q] For the nanowire we employ 
the Thomas-Fermi model (see e.g. Ref. ITT[ ) while for the 
dilute electrolyte we consider a Debye-Hiickel approxi- 
mation (see e.g. Ref. fl2l). We these approximations we 
arrive at the following linear differential equation for the 
induced electrical potential 0, 
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where Atf is the Thomas-Fermi screening length in 
the nanowire domain f2i and Ad is the Debye screen- 
ing length in the electrolyte domains and Q4. The 
electrically insulating oxide-layer domain £^2 is free from 
charges and the potential is simply a solution to the 
Laplace equation. Obviously, the induced potential 
should vanish at infinity and denoting the solution in do- 
main Qi by <f>i and the corresponding dielectric function 
by e, we have the additional boundary conditions 



i = 



re3%i (2b) 



where n is a normal vector to the surface cT^i+i separat- 
ing the neighboring domains fi, and fij+i. Furthermore, 
tr^i+i is the corresponding surface-charge density. In the 
following we consider the case of a surface-charge density 
as accumulated at the surface dfl^^, i.e. oi^ = 02,3 = 
and 03,4 = as- For the dielectric function we have 
£3 = £4. The thickness i of the domain Q3 is physi- 
cally interpreted as the length of the functionalization 
molecules supporting the charge which for £ = resides 
directly on the outside of the oxide-layer of the nanowire. 
In general the oxide layer will be charged unless the an- 
alyte solution equals the isoelectric point of the surface. 
Here we neglect this charge and assume that it may just 
shift the value of the initial carrier concentration in the 
wire. 

In the following we consider the conductance G = Go + 
AG of the nanowire of length L and in particular we focus 



on the conductance change AG due to a finite surface- 
charge density as- In terms of the electron mobility n 
we have G = L^ 1 J n dr en(r) u where n = hq + An 
is the electron density and e is the electron charge. We 
determine the induced charge-density, n = no + An, from 
the Poisson equation, 

-eAn(r) = - ei V 2 <Mr) = —^-<pi(r). (3) 

Formally, the conductance change can now be expressed 
in terms of the induced potential, i.e. 
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where Go = Aeno/j, is the conductance in the absence of 
surface charge. 

Equation ([!]) is our general result for the relative con- 
ductance change expressed in terms of an integral over 
the surface-charge induced potential in the nanowire. 
Combined with Eqs. (|2"a|) and (|2b[) it forms a starting 
point for analytical solutions of simple geometries or nu- 
merical solutions of more complicated geometries, e.g. by 
a finite-element method. 

Before turning to more complicated geometries we first 
consider the case of a long cylindrical nanowire immersed 
in an infinite volume of electrolyte (see Fig. [T]). The 
nanowire has a radius R with an oxide-layer of thickness 
S and due to the cylinder symmetry we may solve the 
problem analytically. For simplicity we focus on the, for 
experiments, highly relevant limit 6 <C R of an ultra-thin 
oxide layer. Furthermore, we consider £ = so that the 
charge resides on the outside of the nanowire oxide layer. 
Solving the linear problem we then arrive at 
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which is easily verified to be a function ranging from zero 
in the limit Ad Atf to unity in the limit Ad 2> Atf- 
We note that in the static limit, water is highly polariz- 
able and £3 ~ 78eo while for silicon ei ~ 12eo. Thus, for 
a fixed no the ultimate sensitivity, Eq. ([T]) with r ~ 1, 
requires a very short Thomas-Fermi screening length, 
but at the same time one would like to benefit from 
the 1/no dependence originating from Go oc no. We 
get Ad = \/ 'e^ksT /2(Ze) 2 co in the range of Ad ~ 1 
- 100 nm for at room temperature and ionic concen- 
trations, Co, for typical physiological electrolytes. For 
example JTo| a realistic salt concentration (Na + Cl~) of 
150 mM yields a screening length of Ad ~ 1 nm thus 
calling for nanowires with densities supporting screen- 
ing at the true nanometer scale. A simple estimate of 

A T f = \J fi 2 £i7r 4 / 3 /m*e 2 n^ 3 also yields Atf ~ 1 nm for 
a carrier concentration of 10 18 cm~ 3 . Taylor expanding 
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FIG. 2: (color online) (a) Sensitivity factor V versus Debye 
screening length Ad- The lower trace: Cylindrical nanowirc 
with radius R and Xtf/R = 0.02. The data points are the 
result of finite-element simulations with a finite oxide layer 
thickness S/R — 0.005. The dashed line shows Eq. ((5}. Upper 
trace: Finite element results for the nanowire in (b) shown 
below, (b) Induced charge-carrier density in a nanowire with 
trapezoidal cross section and Axp/f = 0.02. Aspect ratio: 
w/t = 2, etching defined angle 9 ~ 54, 7°, oxide layer 5/t — 
0.02. Contours show the equipotential lines. 



Eq. (0 in the dilute carrier limit, Atf 3> R, we get 
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so that we arrive at the n 2 ^ 3 scaling in Eq. ([TJ) . 

Figure [2^ illustrates the dependence of the sensitivity 
on the Debye screening in the electrolyte. The lower 
trace is for a nanowire of circular cross section with 
Xtf/R = 0.02. The solid lines shows Eq. © valid for 
a negligible oxide-layer thickness, i.e. S <C R while the 
data points are the result of finite-element simulations 
(Comsol MultiPhysics) taking into account a finite ox- 
ide layer of width 5/R = 0.005 and with €2 = 4. As 
seen, Eq. ((5|) accounts well for the numerical exact re- 
sults. The upper trace shows finite element results for a 
nanowire with a trapezoidal cross section with aspect ra- 



tio w/t — 2, an etching defined angle of 9 ~ 54, 7° corre- 
sponding to the fabricated structure in Ref.d, see Fig. [2b- 
The oxide layer has a thickness 5/t = 0.02 and for the 
Thomas-Fermi screening we have Atf A = 0-02 which 
is somewhat stronger than for the circular case shown 
in the lower trace. Note how the two curves have the 
same overall shape and dependence on the Debye screen- 
ing length, though the stronger Thomas-Fermi screen- 
ing for the upper case makes T approach unity for more 
moderate Debye screening lengths than in the lower case. 
Fig. [2Jd illustrates a typical distribution of the induced 
charge-carrier density in the nanowire and the superim- 
posed contours show the equipotential lines. As expected 
the excess carrier density is induced near the surface of 
the wire supporting the surface-charge density 175. 

Finally, let us discuss the prospects for sensing of point- 
like charges located at a distance I further away from the 
conductor. Obviously, the additional Debye screening in 
the layer of thickness I (see Fig. [T|) will further reduce the 
induced carrier density in the nanowire and in a simple 
picture (neglecting curvature) we would qualitatively ex- 
pect a reduction proportional to exp(— i/Xo)- In the fol- 
lowing we let TV denote the average number of molecules 
absorbed on the wire of length L and we imagine that the 
chain-like functionalization molecule supports a charged 
group, with charge Q, situated at a distance I from the 
surface. Smearing out these charges results in an equiv- 
alent surface-charge density erg = NQ /[2tt(R + £)L] at 
r = R + l. Solving the problem in Eqs. (f2"aj) and pb| for 
a finite i we get T — ► Tg xT with 
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Here, we have used the large-argument exponential 
asymptotes for the Bessel functions. The extra fac- 
tor < < 1 illustrates the additional, close-to- 
exponential, suppression by Debye screening when the 
charge is supported by a functionalization molecule of 
length I. We note that in principle the entering the 
expression for T{ could differ from the Debye screening 
length of the electrolyte e.g. due to the surface function- 
alization. 

In conclusion, we have used Thomas-Fermi and 
Debye-Huckel theory to formulate a simple screen- 
ing model for surface-charge sensing with conducting 
nanowircs. The two screening mechanisms are playing 
in concert and our model illustrates the non-trivial inter- 
play between the nanowire radius R, the Thomas-Fermi 
Screening length Atf, the Debye screening length Ad, 
and the length I of the functionalization molecules. 
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